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You should attempt ALL the questions from Section A, and NOT MORE 
THAN FOUR questions from Section B. Section A carries about 40% of, 
the marks. > 


You may answer questions in any order, writing your answers in the 
answer books provided. Use a separate answer book for cach section. 
At the end of the examination, remember to write your name, student 
number and examination number on the answer books - failure to do so 
will mean that your papers cannot be identified. 
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SECTION A 
(about 40 marks) 


Question 1 


Question 2 


Question 3 


Question 4 


Attempt ALL the questions in this section. 

Not all the questions in this section carry equal marks. 

Write your answers in onc of the answer books provided. Do NOT use the 
same answer book for this section as for Section B. 


Given that x = r cos @ and y = r sin 0, and that u(x, y) is a function on the 
(x, y) plane, express 
ĝu ĝu 
“aty 


K ĝu ĝu 
in terms of r, 0, a and 70° 


A function f with domain R satisfies 
s) = x? 


for 0 < x < 2, and is odd about x = 0 and even about x = 2. Sketch its 
graph for —6 < x < 6 and write down (a) the period of f and (b) the 
formula giving f(x) for 4 < x < 6. 


Match the following three differential equations to the most appropriate 
boundary conditions and domains from the three alternatives supplied, 
so as to give three properly posed problems. 


A dux,y)  , aux, y) 
Tae toa =O 
B Suly) uwy) _ 9 
Ox? ay? $ 
C u(x, y)  „ dulx, y) 
r aye ns 
la ET , 
a Ou _ O<x<1, <»), 
Fz») = 0 
JO =» eye? oe 
b u ĝu pe (0<x<l, <y), 
ax (ls ») = 0455 (0) 1 
u0, y) = y, u(x,0) = 0 
c SY (1,9) = oux O<x<l, 0<y< 1). 


The Fourier series 


i5 r (5) 
= = sin | ~ | cos nax 
n n 3 

n=l 


converges (pointwise) to the sum } when —$ < x < $ and to the sum 
~¢ when 3 < x < 1. To what sum does it converge (a) when x = 4, 


(b) when x = —4? Give the reason for your answer. 


Considered as a serics of functions with domain [— 1, 1], does the series 
converge uniformly? Give the reason for your answer. 
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Question 5 


Question 6 


Question 7 


(a) 


(b) 
(©) 


Question 8 


What is meant by convergence of a finite difference approximation toa 


differential equation? 


A function U(x, t), which satisfies the one-dimensional heat equation 


au _ au 


7 ant? is to be approximated over the mesh x, = th, t} = jk with 


k = 4h", using the following finite difference approximation to the heat 


equation: 
tiger = Wa, + tar] - 
Given that the true solution U satisfies 


Ui j+i = HO. + Gi41.4) + O(k?) 


where U, ; means U(x,,t)), deduce that the approximation method is 


convergent. 


Laplace's equation in two-dimensional polar coordinates (r, 8) is 
ðu lau 1 au f 
at ta t ae 

The method of separation of variables yields as a solution 
u = (Ar* + Br-*C cos k0 + D sin k0) 


=0 (r>0,0 <0 < 2x). 


where A, B, C, D, k are constants and k > 0. For a solution of this type 
what conditions on these constants are implied (a) by the condition that u 
must be single-valued, and bounded in the neighbourhood of the origin, 
(b) by the condition that u must be single-valued and bounded for 


arbitrarily large r? 


This question is about the system of equations 
10x, — x, =1 
=x; + 10x- x3 =2 n3 

— xXx + 10x; =3. 


Starting with the approximation x, = x, = x, = 0, calculate the next 


approximation using the Jacobi method. 

Same as (a), but using the Gauss-Seidel method. 

Given that the eigenvalues of the Jacobi iteration matrix 
0 wy 0 
wo > 
0 %0 


arc 0, ys V2, — ty V2, give an estimate of the factor by which the accuracy 
of the approximate solution is increased at each application of the Gauss- 


Seidel method. 


By comparing with a suitable constant-coefficient equation, obtain a 


Positive lower bound for the lowest cigenvalue A of the problem 
2 
TD — ax) + dia) = 0 (l <x <2) 


ul) = (2) =0. 
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SECTION B 
(about 60 marks) 


Question 9 (i) 


(ii) 


(iii) 


(iv) 


Question 10 (i) 


ii) 


(ii) 
(iv) 
(vy) 


Attempt at most FOUR questions from this part. 
Write your answers in one of the answer books provided. Do NOT use the 
same answer book for this section as for Section A. 


Explain what is meant by the word characteristics in connection with 
hyperbolic equations. 
Find the characteristics of the equation 
ðu ĝu 3u 
det — 9) Seay + 23 = = 0 (x, yeER). 


Hence, or otherwise, find a change of independent variables which reduces 


9 


Ch is 
this system to the standard form a = 0 (it is not necessary to carry 
out the actual reduction). 
a 
If the values of u and x are specified along the straight line 


y=mx+c, 


for what values of the constants m and c would it be impossible to find a 
unique solution? Give the reason for your answer. 


Classify the equation 


2, 2 
dul, ) qy Pula 1) 4x2 THD _ 9 (x, teR) 


Ox? Ox ôt Oe 
as elliptic, hyperbolic, or parabolic. 


How many characteristics of this equation pass through a typical point in 
the (x, t) plane? 


Find the characteristics. 
Find a coordinate system in which the equation takes its standard form. 
> 


Carry out the transformation to standard form. 
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Question’ 11 


Question 12 


The diagram shows the cross-section of a pipe which is uniform in the 
direction perpendicular to the paper. The inner contour is rectangular, 
Write down a differential equation and boundary conditions-which model 
the steady flow of a viscous fluid along this pipe, defining any new symbols 
introduced (apart from those in the diagram) and Stating what assumptions 
you are making about the direction and about the space and time depen- 
dence of the velocity vector. 


Use the second extremum principle for viscous flow, with a trial vector 
field of the form v* = (ax, By) where a, B are constants, to obtain an upper 
bound on the volume of fluid per unit time delivered by the pipe. What 
values of a, $ give the best such bound? 


Let vo” denote the trial vector field (not necessarily of the form (ax, By)) 
which gives the best possible upper bound in the second extremum 
principle. Express the components of vọ* in terms of the velocity of the 
fluid in the pipe and its derivatives. 


Describe one explicit and one implicit finite-difference scheme for solving 
the problem 


U(x, t) _ U(x, t) U(x, t) 
Ta Og t r a 


U(0,1t)= 0 
aU.) o (2 0) 
Ox 
U(x,0) =x? -2x (<x). 
Compare and contrast your two schemes with respect to stability, accuracy, 


and ease of computation, including detailed reasons for any assertions, 
you make. i 


Sera Turn over 


Question 13 (i) 


(i) 
(iii) 
(iv) 


Question 14 


Gi) 


(iii) 


Bring the equation 


2 k 
Fe + u Ox <n 


to self-adjoint form. 


Find the Green's function G(x, č) for the problem of solving the above 
equation with the boundary conditions 


4(0) = u(1) = 0. 
Verify, if you have not already done so as part of the calculation of 
G(x, ¢), that G(x, ) is continuous and has the correct discontinuity in its 
derivative at x = &. 


Indicate briefly, and without doing any calculations, how you could use 
this Green’s function to solve the given problem. 


Let D denote the semicircular region of the (x, y) plane shown in the 
diagram. If u is a function defined on D, satisfying 5 = Owhen y = 0, 


show that the corresponding function for polar Coordinates (which we also 
denote by u) satisfies ni regy 


Ou(r, 0) = Our, n) 0 
; -a nr a 
Reformulate the following problem in Polar coordinates: 


ðu au 
oe t Rt F (Œœ, y) € D) 


(O<r<1). 


u=0 on the boundary curve ABC (1) 


= = 0 on the boundary AOC 


where F is an arbitrary continuous function on D. 


The Green's function for the problem (1), expressed in polar coordinates, 
has the form 


: r? — 2rp cos (0 — $) + p? 
G(r, 8; p, $) = ain fiz = 2rp cos (0 — $) + i} 


r? — 2rp cos (0 + ore} . 


tong P7 — 2rp cos (0 + $) +1 


Evaluate the constants a and b. 
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Question 15 


6) 
Gi) 


(iii) 


Question 16 


Question 17 


@ 
Gi) 


A Dirichlet problem for the equation 


U(x, y) | U(x, y) _ U(x, y) 
ax + dy? + Fx =0 


is to be solved approximately, using a square mesh of side h and the 
finite-difference scheme 

guig — (1 + tdi — (L te — tga — uj- =O 
where u, j is the approximation to U(x, yj and x, = xo + ih, Yy = Yo + jh. 


For what range of values of h does the resulting system of simultaneous 
equations satisfy a maximum principle? 


Given that A satisfies the condition obtained in (i), use the maximum 
principle to prove that the system of equations has a unique solution. 


Prove that the finite difference scheme is a consistent approximation to the 
differential equation. 


Use the method of separation of variables to obtain the general solution 
in series for the equation 


ðu lôu du 
Ttr t aA 
where u = u(r, z, t) and a, b, c are positive constants, subject to the 
boundary conditions: 
u(r, z, t) bounded in the neighbourhood of r = 0, 
ur, z,t) = Oifr =a, or ifz =O0orz = b, or if t = 0. 
Express your answer in terms of the zero-order Bessel function Jo and its 
zeros j,°. 


=357 O<r<a0<z<6,1>0) 


It is proposed to use the finite-difference scheme given by the 
approximations 

Ou 1 g 
E =p 67m, > 


Chin 1 
“Zyl & a lað, uge + (1 — 2a) 5,2u,) + adaugi] 5 
where a is a constant, to calculate an approximate solution of 


a U(x, 


over the mesh x, = Xo + ih, t) = tọ + jk, where u,, stands for the 
approximation to U(x; t)). 
Write out the finite difference scheme. 


What does the von Neumann method indicate about the stability of this 
scheme when a = $? 
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